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The recent discovery of topological insulators has revived interest in the band topology of insu- 
lators. In this work, we extend the topological classification of band structures to include certain 
crystal point group symmetry. We find a class of three-dimensional "topological crystalline insula- 
tors" which have metallic surface states with quadratic band degeneracy on high symmetry crystal 
surfaces. These topological crystalline insulators are the counterpart of topological insulators in 
materials without spin-orbit coupling. Their band structures are characterized by new topological 
invariants. We hope this work will enlarge the family of topological phases in band insulators and 
stimulate the search for them in real materials. 

PACS numbers: 73.20.-r, 73.43.-f 



The topology of band structures is important in the 
study of topological phases of matter. Historically the 
quantization of Hall conductance in the integer quan- 
tum Hall effect was explained by the Thouless-Kohmoto- 
Nightingale-den Nijs integer (also known as Chern num- 
ber) of occupied energy bandspQ. Recently the topolog- 
ical classification of spin-orbital coupled band structures 
with time reversal symmetry has played a key role in the 
theoretical prediction of topological insulators |2j|4]. The 
subsequent development of topological band theory [5], 
combined with realistic band structure calculations, has 
proven useful in the material search for topological insu- 
lators. 

Inspired by the discovery of topological insulators, 
the classification of band structures has been extended 
to other discrete symmetry classes such as particle-hole 
symmetry [6 8J, which leads to a rich family of topo- 
logical phases such as topological superconductors E]. 
More recently, the classification of magnetic insulators 
with certain magnetic translation symmetry has been 
studied[ini ITT] . Finding these phases in real materials 
is interesting and challenging. 

The purpose of this work is to extend the classifica- 
tion of band structures in a different direction — to in- 
clude crystal point group symmetries. We introduce the 
notion of "topological crystalline insulators" , which can- 
not be smoothly connected to a trivial atomic insulator 
when time reversal (T) symmetry and certain point group 
symmetry are respected. Unlike time reversal symmetry, 
crystal symmetries can be broken by sample surfaces. As 
a result, a low-symmetry surface of a topological crys- 
talline insulator does not have robust surface states. This 
motivates us to study a class of three-dimensional topo- 
logical crystalline insulators which have four-fold (C4) or 
six-fold {Cq) rotational symmetry. As we will show, its 
(001) surface, which preserves the rotational symmetry, 
supports gapless surface states. These topological crys- 
talline insulators are the counterpart of topological in- 
sulators in materials without spin-orbit coupling. Instead 
electron's orbital degrees of freedom play a role similar to 
spin. Unlike the linearly dispersing Dirac surface states 
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FIG. 1: (a) tetragonal lattice with two atoms A and B along 
the c-axis in the unit cell; (b) the Brillouin zone and four high 
symmetry points. 



of topological insulators, the (001) surface states of topo- 
logical crystalline insulators have quadratic band degen- 
eracy protected by time reversal and discrete rotational 
symmetry [TH [T3] . 

The outline of this paper is as follows. First we intro- 
duce a simple tight-binding model in a tetragonal crystal 
with C4 symmetry. We explicitly show that gapless sur- 
face states exist on the (001) crystal face. The topological 
stability of surface states suggests a nontrivial phase in 
this model. Next we define a new Z2 topological invariant 
for generic time reversal invariant band structures with 
C4 or Cg symmetry in three dimensions. This estab- 
lishes the topological crystalline insulator phase. Finally 
we discuss the electronic properties of quadratic surface 
bands. 

Tight-Binding Model: Consider a tetragonal lat- 
tice with a unit cell of two inequivalent atoms A and B 
along c-axis, as shown in Fig. la. The 3D crystal can 
be viewed as a stack of bilayer square lattices in the ah 
plane. We now introduce a tight-binding model to de- 
scribe the band structure of electron's p orbitals (or d 
orbitals, see the next paragraph). In particular, we are 
interested in the energy bands derived from p x and p y or- 
bitals. We assume that these bands do not overlap with 
the p z bands, and construct a tight-binding model from 
Wannier functions with the same symmetry as p x and p y 
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orbitals. The Hamiltonian H consists of intra-layer hop- 
ping H A and H B , as well as inter- layer hopping H AB : 

H = / ] H A + H B + H„ B , 

n 

Hn B = ^(T i -r j )[^c] la {T i ,n)c Ba (T j ,n) + h.c.] 

i,j a 

+ t' z Y.Y} c ^ a (^ n )cB a {v l ,n + l) + h.c] (1) 



Here each site is specified by (n,r, a): n labels the bi- 
layer unit cell along c-axis; r = (x,y) labels the a6-plane 
coordinate; a = A, B labels the sublattice. a, (3 label the 
p x ,Py orbital. Two types of p-orbital hopping terms ap- 
pear in H. The intra-layer hopping in H a is of c-bonding 
type: it depends on the relative orientation of p-orbital 
and the hopping direction e* 3 ' = (r^ — ij)/|rj — r j\- The 
inter-layer hopping in H AB is orbital- independent. Writ- 
ten in this form, the tight-binding Hamiltonian mani- 
festly preserves crystal symmetries. The hopping ampli- 
tudes t a , t' and t' z between two sites depend on their 
a6-plane distance — Vj. Throughout this work, we as- 
sume spin-orbit coupling is negligible, so that electron's 
spin index is omitted. 

We emphasize that the form of the Hamiltonian H is 
entirely determined by crystal symmetry. Because d XZjVZ 
orbitals transform in the same way as p XiV orbitals un- 
der O4, ([!]) also applies to materials with these d-orbital. 
Therefore ([I]) is potentially revelant to a large class of 
materials including transition metal compounds with t2 g 
orbitals near Fermi energy. 

To obtain a minimal model for topological crystalline 
insulators, we include the nearest and next-nearest neigh- 
bor intra-layer hoppings in H a with the amplitude t\ and 
£§) as well as nearest and next-nearest neighbor inter- 
layer hoppings in H AB with the amplitude t[ and t' 2 . 
The corresponding Bloch Hamiltonian -ff(k) is obtained 
by Fourier transform: 



H a (k) 



= 2tl 



H A (k) 
H AB ^(k) 

cos k x 




H AB (k) 
H B (k) 

V 

cos ky 



2*S 



H AB (k) = [t\ 



cos k x cos k y sin k x sin k y 
sin k x sin k y cos k x cos k y 

- 2t' 2 (cos k x + cos k y ) 



]!■ (2) 



The band structure is shown in Fig. 2a for the following 
set of parameters: tf = -tf = l,t A = —t B = 0.5,^ = 



2.5,^ = 0.5, t' 



2. We have checked that the energy 



gap is finite everywhere in the Brillouin zone. As long as 
the energy gap does not close, the system remains in the 
same topological class within a finite parameter range. 
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FIG. 2: (a) bulk band structure of the tight binding model 
along high symmetry lines; (b) surface states with quadratic 
band touching exist on (001) face. The tight-binding param- 
eters are shown in the text below Eq.|2]|. 



To study surface states, we solve H in a slab geome- 
try. We find that the existence of gapless surface states 
crucially depend on the surface termination. Consider 
the high symmetry (001) surface, which preserves the C4 
symmetry. Here surface states exist and transverse the 
whole energy gap as shown in Fig. 2, leading to a 2D 
surface metal. Note that surface states are doubly de- 
generate at M = (7r, 7r): one in p x orbital and the other 
in p y orbital. Because M is a fixed point under four-fold 
rotation, the doublet form a two-dimensional irreducible 
real representation of C4 symmetry, as a result of time 
reversal symmetry for electrons without spin-orbit cou- 
pling (i.e., spinless fermions). Surface band dispersion 
near M can be obtained from k ■ p theory. We introduce 
a pseudospin a z = ±1 to label the p x and p y orbital of 
the doublet. In this basis, C4 rotation is represented by 
e «Tj,7r/4 anc j tj me reversal operator T is represented by 
complex conjugation. Up to a gauge transformation, the 
form of the k ■ p Hamiltonian % is dictated by symmetry: 



l~L{k x , ky 



k 2 

Ky 



2mn 



2mi 



k x k y 
2m 2 



(3) 



Perturbations which break either C4 or T symmetry can 
open up an energy gap and destroy the protected sur- 
face states. This can be seen explicitly by adding a C4- 
breaking term M\k x o y+Mio z or a T-breaking term Ma y 
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to %. Similar quadratic bands have been recently stud- 
ied in 2D photonic crystals [12] and fermion models [13]. 
Our tight-binding model does not have gapless surface 
states for other surface terminations which break the C4 
symmetry. 

When T and C4 symmetry are preserved, the surface 
states shown above are topologically protected. This can 
be understood by considering the surface band dispersion 
along f M. Within our model, surface states (if present) 
must be doubly degenerate at f and M. As a conse- 
quence, there exists two distinct types of surface band 
connectivity along TM, having an even or odd number 
of band crossings at Fermi energy respectively. Surface 
states with even crossings are fragile: they can be pushed 
out of the energy gap by changing the surface potential. 
Surface states with odd crossings have a Fermi surface 
which encloses T an odd number of times. They can- 
not be removed. These two types of surface states imply 
the existence of two topologically distinct phases of time- 
reversal-invariant band insulators with four-fold symme- 
try. The reasoning here closely follows the previous study 
of topological insulators 14]. 

Z2 Topological Invariant: We now show that a Z2 
topological invariant Vq = 0, 1 characterizes the band 
structure of 3D time-reversal-invariant insulators with 
four-fold rotation symmetry (without spin-orbit cou- 
pling) .^o = corresponds to a trivial phase adiabatically 
connected to an atomic insulator, isq = 1 corresponds to 
a topological crystalline insulator which has gapless sur- 
face states on the (001) surface. 

We start by examining the symmetry property of Bloch 
wavefunctions of occupied bands: |^„(k)) = e ik ' r |u„(k)), 
where |u„(k)) is the cell-periodic eigenstates of the Bloch 
Hamiltonian H(k) = e~ ik r i/e ik r . Here the unit cell is 
chosen to be invariant under C4 rotation around the z 
axis. H(k) satisfies 



H{k x ,k y ,k z ) — UH(ky, k x ,k z )U 
H{k x ,ky,k z ) = TH(—k x , ky, k z )T 



(4) 



Here the unitary operator U — e lLz7r ^ 2 implements C4 
rotation within the unit cell (L z — xP y — yP x is an- 
gular momentum operator). The anti-unitary operator 
T = K (complex conjugation) implements time rever- 
sal transformation for spinless fermions, with the prop- 
erty T = 1. As a result, time reversal symmetry by 
itself does not guarantee a two-fold degeneracy. How- 
ever, the combination of time reversal and four-fold ro- 
tational symmetry can lead to protected degeneracies at 
four special momenta T = (0, 0,0), Af = (tt,tt,0),A = 
(tt,tt,tt),Z — (0,0, 7r) in the 3D Brillouin zone (Fig. lb). 
At such a high symmetry point k;, H (kj) commutes with 
U, so that the energy bands \u n (ki)) are eigenstates of 
four-fold rotation with possible eigenvalues 1,-1, i and 
—i. Moreover, because £f(kj) is real, energy bands at 
ki with ±i eigenvalues are guaranteed to be degenerate, 



forming a two-dimensional irreducible real representation 
of C4. For example, such doublet bands can derive from 
(PxiPy) orbitals or (d xz ,d yz ) orbitals. From now on, we 
consider a set of energy bands |u n (k)), n = 1, ...2N which 
are doubly degenerate £? 2 n-i(kj) = E 2n 0<-i) at r, M, A 
and Z. As we will see below, only these doublet bands 
admit a Z2 classification. 

The Z2 topological invariant v is defined in terms of 
the electron wavefunctions: 



(-ir 

(— l)^i k a 
A k 



(_1)^M(_ 1 )«A Z) (5) 

, /" k2 „ „ N PfMk 2 )] 
-i£>i(k)|ac|«i(k)>, 



w mn (ki) = (u m (k l )|{7T|M„(k i )). 

„4k is the U(l) Berry connection. u;(kj) is an anti- 
symmetric U(2N) matrix because [H{ki),UT] = and 
(UT) 2 = —1. Pf stands for the Pfafnan. We now specify 
the integration path in ^ . For u^m , the integral is along 
an arbitrary line connecting T and M which lies within 
the 2D plane k z — in the Brillouin zone as shown in 
Fig. lb. Likewise, the integration path for vaz li es m the 
plane k z = ir. 

Both A and w depend on the choice of basis \uj(k)), 
but we now prove that ^kik 2 is gauge invariant. Two 
different basis |u„(k)) and |u„(k)) are related by a U (2N) 
gauge transformation: 



,(k)>=e„ m (k)|u TO (k)) ) GeU(2N) 



(7) 



This leads to the corresponding gauge transformation of 
A and w: 



A = A- iTr^-^a] =A-tdk log(det[£]) 



(8) 



Using the identity Pf[A T A/A] = det[A]Pf[M], it is 
straightforward to verify that ^kik 2 — is gauge in- 

variant. 

We now further show that (— l) Uk i k 2 = ±1 is a Z 2 
quantity. It follows from Q that at the two 2D planes 
k z = and k z = ir, H(k) has the symmetry property: 
H(k) = EH(k)E, S = U 2 T. This allows us to choose a 
real gauge along the integration path to evaluate Q : 



3|w ro (k)) = -|w m (k)), 



(9) 



Because H is anti-unitary, A — everywhere along the 
integration path. So in this gauge (— l) 1/ [ kl ' k2 l reduces to 



(_!)^ lk2 =Pf[ u ,(k 2 )]/Pf[ u ,(k 1 )]. 



(10) 



Because |u m (kj)) belong to the two-dimensional repre- 
sentation of C4 group, we have [/ 2 |it m (k;)) = — \u m (ki)), 
so that the gauge condition (|9| at k^ is equivalent to 
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T|u m (kj)) = \u m (ki)), i.e., the wavefunction is real. Un- 
der this reality condition, the matrix iu(kj) simplifies to 

«w(kj) = (M m (k i )|C/|u„(k i )). (11) 

Now by choosing a particular kind of real basis 
|w2m(kj)) = [7|tt2 TO -i(k,)), u?(kj) reduces to a standard 
form w° = e © e... © e, which is a direct sum of iV 2 x 2 
Levi-Civita tensors. This means that in a generic real 
basis, w(ki) can be written as 

wQ*) = O T (k i )w O(k l ), 0{ki) e 0(2N). (12) 

So we have 

Pf[to(ki)] - det[O(k i )]Pf[w ] - det[0(k i )] - ±1. (13) 

This proves that (— l) lyk i k 2 = ±1 is a gauge invariant Zi 
quantity. So both v^m and z^z are Z 2 topological in- 
variants which characterize the band structures in the 2D 
momentum space k z = and k z — 7r respectively. Their 
product is defined for 3D time-reversal-invariant in- 
sulators with C4 symmetry, and determines the surface 
state property: gapless surface states exist on (001) face 
only when f = 1. This is verified in the tight-binding 
model ([I]). The relation between i>o and i^tm^az is anal- 
ogous to that of strong and weak Zi index in 3D topo- 
logical insulators |15l4T7] . 

Generalization to Crystals with Cq Symmetry: 
The study of topological crystalline insulators with C4 
symmetry can be generalized to hexagonal crystal struc- 
tures with Cq symmetry. In that case, (001) surface 
states have quadratic degeneracy points either at f or 
two equivalent K. The topological invariant is similarly 
defined by Eq. |6]) from the electron wavefunctions on the 
lines TK and AH in the Brillouin zone. A concrete tight- 
binding model can be obtained by placing ([I]) in a layered 
triangular lattice. The details will be reported elsewhere. 

The fact that the above Z2 invariant is only defined 
for doublet bands illustrates the important interplay be- 
tween symmetry representation and topology of energy 
bands in solids, a feature absent in other known classes 
of topological insulators. We are not aware of a system- 
atic mathematical approach to classifying vector bundles 
(e.g., electron wavefunctions over the Brillouin zone) with 
given representations of group action (e.g., C4 rotation) 
at fixed points (e.g. high symmetry momenta kj). We 
hope the concrete topological invariant ^ can stimulate 
continuing research on topological crystalline insulators 
with other lattice structures. 

The occupied bands of real materials usually have 
both doubly degenerate and nondegenerate singlet com- 
ponents. Provided that the doublet bands are energeti- 
cally separable from the singlet ones (i.e., no band cross- 
ing between them), the Z2 invariant remains well-defined 
for the former. The coexistence with singlet bands can in 
principle weaken the stability of surface states (see sup- 
plementary material), although this senario seems un- 
likely to happen in reality. 



Surface States with Quadratic Degeneracy: 

Quadratic degeneracies in 2D band structure have at- 
tracted considerable interest recently. The k ■ p Hamil- 
tonian ([3]) is widely used as a low energy approxima- 
tion for the band structure of bilayer graphene near ±K 
points. However, unlike topological crystalline insula- 
tors, the quadratic degeneracy there is not protected by 
symmetry. Instead, the band dispersion very close to ±K 
becomes linear due to trigonal warping effects. 

Electron interactions can drive 2D bands with 
quadratic degeneracy into a variety of interesting bro- 
ken symmetry phases including quantum anomalous Hall 
state, nematic phase and etc|13j. The competition be- 
tween different ordered phases is of great interest in bi- 
layer graphene. It remains to see what happens in the 
single-valley quadratic surface states of topological crys- 
talline insulators. 

It would be interesting to generalize the concept of 
topological crystalline insulators to interacting systems. 
For example, crystal symmetry is known to play an 
important role in the topological classification of spin 
liquids [15], spin chains [H?]. and Mott insulators[5D]. A 
unifying theory of symmetry protected topological phases 
remains to be developed. We also note that analogs of 
topological crystalline insulators can be potentially real- 
ized in photonic crystals. 

Part of this work was done in early 2009 at University 
of Pennsylvania under the support of NSF grant DMR- 
0906175. We thank Charlie Kane and Jeffrey Teo for 
very helpful discussions, and acknowledge support from 
the Harvard Society of Fellows. 
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SUPPLEMENTARY MATERIAL 



W = 




We study the stability of surface states with quadratic 
degeneracy when the doublet bands coexist with singlet 
band in the bulk. The latter may come from, e.g., p z 
orbitals so far neglected in our tight-binding model 
on tetragonal lattice. Provided there is no band crossing 
between doublet and singlet bands, the invariant v 
for the former remains well-defined. However, as we now 
demonstrate, the quadratic surface states of a nontrivial 
v = 1 phase can in principle be damaged by hybridizing 
with additional singlet surface band. 

We assume that the singlet band in the bulk gives rise 
to a branch of unprotected nondegenerate surface state 
with the usual dispersion E 3 (k) = e 3 + fc 2 /(2m 3 ). e 3 
is the energy difference between the doublet and singlet 
surface state at k = 0. This extra surface state can hy- 
bridize with the original surface states with quadratic 



Here % is the 2x2 Hamiltonian defined in Eq.Q. A 
is the coupling strength between them responsible for 
hybridization. For 637713 < and A = 0, singlet and 
doublet surface bands cross at k ^ 0. A finite coupling 
A 7^ turns the band crossing into an avoided crossing 
and opens up an energy gap. The resulting surface states 
cease to be protected: they can be smoothly pushed into 
the bulk continuum. This shows that the intervention 
from singlet bands can in principle destroy the gapless 
surface state with quadratic degeneracy. However, the 
above scenario requires the coexistence of singlet surface 
state under particular conditions, which seems very non- 
generic in real materials. Therefore we are optimistic 
that surface states with quadratic degeneracy on topo- 
logical crystalline insulators are robust for practical pur- 
poses. 



